Transport parameters in neutron stars from in-medium NN cross sections 
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We present a numerical study of shear viscosity and thermal conductivity of symmetric nuclear 
matter, pure neutron matter and /3-stable nuclear matter, in the framework of the Brueckner the- 
ory. The calculation of in-medium cross sections and nucleon effective masses is performed with a 
consistent two and three body interaction. The investigation covers a wide baryon density range as 
requested in the applications to neutron stars. The results for the transport coefficients in /3-stable 
nuclear matter are used to make preliminary predictions on the damping time scales of non radial 
modes in neutron stars. 

PACS numbers: 97.60. Jd,21.65.-f,26.60.-c, 



I. INTRODUCTION 



Neutron stars are a unique laboratory for studying the equation of state of nuclear matter at high density and 
isospin asymmetry, beyond the ranges typical of heavy ion collisions and nuclei far from stability. But the observation 
data are still far from uniquely constraining the properties of exotic states of nuclear matter. The detection of 
gravitational waves could push forward the neutron star physics. As the gravitational waves drive the instability of 
I neutron star oscillations, including non radial modes, the possible damping mechanisms have to be investigated to 
justify the existence of rapidly rotating stars . Good candidates are the viscosity and the thermal conductivity of the 
Q neutron star constituents. The important role that these parameters can play is the main reason for the uninterrupted 
i—{ •■ interest in the transport theory of dense matter over the last three decades Microscopic models of nuclear matter, 
CZ3 , i.e. based on bare interactions, have been faced with the interpretation of the equilibrium properties of neutron stars, 
, such as their structure [3 and the onset of superfluidity [1]. The same theoretical models have been also applied to 

[ calculate effective mass^JJand medium renormalized nucleon- nucleon (NN) cross sections [5*1, both extensively used in 
^-H ■ the transport-model simulations of heavy ion collisions. But only recently these two quantities were redirected to 
^ ' describe the transport properties of nuclear matter. These quantities are in fact the main ingredients to calculate 
the viscosity and thermal conductivity coefficients in neutron stars. The extension of such calculations to regimes of 
- high density and isospin, needed in the study of neutron matter core, together with the interpretation of their inner 
, structure, is a main challenge for the current microscopic theories of the nuclear matter. 
• ' So far all calculations based on microscopic many body approaches with realistic interactions have shown that 
I the medium effects result in a noticeable suppression of the NN cross sections UATAfQ and, as a consequence, in 
an enhancement of both shear viscosity and thermal conductivity. Calculations in dense nuclear matter have been 
performed in the T-matrix approach (6|, in the correlated basis function (CBF) approach[7|, and Brueckner- Hartree- 
Fock (BHF) approximation 0. Calculations of the transport coefficients in /3-equilibrium nuclear matter have been 
performed with the equation of state (EoS) of asymmetric nuclear matter from the variational approach fTl| . and 
from CBF approach 9]. Despite the overall agreement on the medium effects, different approaches could give different 
predictions for transport parameters in /3-equilibrium nuclear matter since they differ from each other in the different 
isospin dependence that becomes more and more visible in higher density range. The results should be insensitive to the 
choice of the two body potential because all realistic potentials used in the microscopic calculations are accommodated 
on the experimental NN scattering phase shifts. But a dependence is expected on the three body force, especially 
at the high density, where its influence on the EoS is dominant. Therefore, it seems interesting to compute the 
shear viscosity and the thermal conductivity of nuclear matter, especially /3-stable nuclear matter, in a wide range of 
densities needed for the study of the neutron star core. This will be done within the BHF approximation. The latter 
embodies, within a unified meson-exchange model framework, two and three body forces. The meson parameters of 
the two-body realistic interaction (the Bonn B [l^ in our case), which fit the experimental NN scattering phase shifts 
in the vacuum, are also adopted to describe the three-body force (l2j . In this paper the numerical results will be 
tested, within the simple model of constant density neutron star, on the calculation of the dissipation time scales to 
be compared with the time scale of emitting gravitational radiations. 
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II. TRANSPORT PARAMETERS 

The transport parameters of Fermi liquids were derived by Abrikosov and Khalatnikov (AK) from the Landau 
kinetic equations for a muhicomponent svstemp3j 

where fi (r , p) is the quasiparticle distribution of the component i , ep the quasiparticle energy, and lik the colhsion 
integral between particles of the components i and k. From the linearization of the kinetic equations, the shear 
viscosity r/ and thermal conductivity k can be extracted in the AK approximation. The exact expressions, obtained 
by Brooker and Sykes after revisiting the resolution of the kinetic equations [TBI [l6j. are written as follows 

VT' - ^pvlW{p)C{X) (2) 
kT = lvj,p^W{p)H{fi), (3) 



where 



1 f^'" .„ Z"^" de 1 



where p is the density, vf = pp/m* is the Fermi velocity and m* the effective mass. C(A) and H{p) are correction 
factors corresponding to the exact solution (see for details [H, [IB])- The key quantity is the in medium cross section 
here expressed in terms of the energy E in the laboratory frame and the scattering angle 9 in the center of mass 
frame. The upper limit of the energy integration in the average cross section is four times the Fermi energy of free 
Fermi gas ep = p'pl'^fn due to the approximation of restricting the nucleon excitations around the Fermi surface. As 
described below, the in medium cross sections are calculated within the Brueckner theory. 

III. IN MEDIUM CROSS SECTIONS FROM THE BRUECKNER THEORY 

In the interior of a neutron star the hadron density can reach values several time the nuclear matter saturation 
density, so that the nucleon-nucleon (NN) collisions are expected to be deeply affected by the surrounding medium 
and the corresponding cross sections can be quite different from those in free space. There are two main medium 
effects: first, the NN scattering amplitude is dominated by the S-wave components of the effective interaction so that 
a flattening of the angular distribution is expected in comparison with the free cross section, which, in the center of 
mass frame, is peaked in the forward and backward directions ; second, the level density in the entrance and exit 
channels gets reduced by the strong medium renormalization of the effective mass. 

Both effects can be well described in the framework of the selfconsistent Brueckner theory. In the last years the 
latter has made a remarkable step forward by means of the three-body force, introduced not only to reproduce the 
empirical saturation properties of nuclear matter but also to extend to high density the calculations. The Brueckner 
theory with two and three-body force is described elsewhere [13, [3 ■ Here we simply give a brief review of the BHF 
approximation, adopted for the present calculations. The starting point is the reaction G-matrix, which satisfies the 
Brueckner-Bethe-Goldstone (BBG) equation, 

G(i:j) = WNN + WNN > G(a;), (5) 

f-f w - Efei - efc2 

where ki = (fci, di, r^), denotes the single-particle (s.p.) momentrmi, the ^-component of spin and isospin, respectively, 
and w is the starting energy. The G-matrix, the Pauli operator Q and the s.p. energies Cfc = fc^/2m -I- Uk depend 
separately on the neutron and proton densities. The interaction wnn given by 

^^NN = ^2'^"'= + Vf, (6) 
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where Vj"^''® is the bare two-body force and V^^ is the three-body force averaged on the third particle as foUows 

{rir2\VfiTM%') = 

^Tr J2 I dradfg'C, (r^) [1 - i^{r[^)] [1 - ii{r'^^)] (7) 

xT4^3(f/f2'r3'|fl?^2r=*3)(/'n(r3)[l - '7(?'13)][1 - ??(?'23)]- 

Since the defect function 77 (r) (1 — 77 (r) is the correlated two body wave function) is directly determined by the 
solution of the BBG equation (17J, V^^ must be calculated self-consistently with the G-matrix and the s.p. potential 
Uk on the basis of the self-consistent BBG equations. It is clear from Eq. Q that the effective force rising from the 
3BF in nuclear medium is density dependent via the defect function. A detailed description and justification of the 
method can be found in Refs. [171 [l8l|. 

In the present calculations the Bonn B potential was adopted as Vj"^''*' [l^. Besides being a realistic interaction 
fitting the experimental the NN scattering phase shifts, it has the advantage of being built up in terms of meson 
exchange so as the three body force. Therefore the choice of the same meson parameters, that is masses, coupling 
constants and cutoffs, provides a unified treatment of two and three-body forces as mentioned before. Details and 
results with this interaction are presented in Ref. [l2| . In FigH^ the equations of state of pure neutron matter (PNM) 
and symmetric nuclear matter (SNM) are plotted. The EoS from the correlated basis theory (CBT) Q and variational 
chain summation (VCS) approach 19] are also plotted for a comparison , that will be useful for the next discussion 
of the transport parameters. A strong deviation from the Brueckner calculation can be easily observed at increasing 
density, which is much more sizeable for the symmetry energy plotted in FigUlD. The symmetry energy affects the 
neutron star transport properties, because it essentially determines the isospin composition of the core. 




FIG. 1: (color online). (a) Energy per particle in pure neutron matter and symmetric nuclear matter. (b) Symmetry energy, 
calculated as the difference between the neutron the symmetric nuclear matter energies. The diamonds represent the results 
obtained using the correlated-base function (CBF) approximation with a Fermi gas states 01, the dashed lines represent the 
results obtained using the variational chain summation (VCS)[l9[. 



In the Brueckner theory the in-medium NN cross section is obtained replacing the T- matrix with the G-matrix and 
the in vacuum level density with the in-medium one. This definition is supported by the property of the G-matrix to 
go over into the T-matrix in the zero density limit, as it can be seen from BBG equation ( Eq.([5])). 

In the case of pure neutron matter the neutron-neutron cross section in the center-of-mass frame is written 

'^™(i^,^) = Y£i E \Gls'SO) + {-lfGls'M-e)\^- (8) 

The prefactor is the square of the level density at the Fermi energy, where the AK approximation assumes the particle 
transitions to mainly occur. The medium renormalization of the nucleon mass to values much less than the unit. 
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FIG. 2: (color online). Density dependence of effective mass mj^/m in pure neutron matter (blue solid line), symmetric nuclear 
matter (black dashed line) and the /3-stable nuclear matter (pink dotted line), respectively. 



as shown in Fig [21 reduces the level density with respect to the free Fermi gas value. As a consequence also the in- 
medium cross section turns out to get reduced. The additional medium effect, which is incorporated in the G-matrix, 
is due to the Pauli blocking which prevents the particles to scatter into occupied states. In Fig. [3] (panel (a)) the 
in- medium neutron- neutron cross section at the laboratory energy E ~ 100 MeV, obtained from Eq.([8|), is compared 
to the corresponding free one. As mentioned above, the differential cross section becomes more and more isotropic 
at increasing density. In addition, its magnitude is reduced, that is a common feature to all predictions based on 
microscopic approaches. In all curves the angular distribution is symmetric around 90 degrees, since the cross section 
section is antisymmetrized for identical particles. In the panel (b) the neutron-neutron cross section in various baryon 
environments is reported. It is worth noticing that it is completely isotropic in /3-stable matter. 

In nuclear matter, besides the neutron-neutron scattering, neutron-proton scattering must be considered. In this 
case the cross section is 

-"p(^'^) = l£^ E IG'ls'M'- (9) 

In FiglSj; the corresponding cross sections are depicted for symmetric nuclear matter and /3-stable matter. A common 
feature to all cases is the medium suppression, but a difference is to be noticed in the angular distribution. The 
enhancement in the backward direction of cr„p is a signature of the anisotropic behavior of the tensor force in the SD 
channel, giving the dominant contribution to the neutron-proton interaction [20l [2l| . 

In the core of neutron stars the neutron and proton composition is determined by the condition of equilibrium with 
leptons (electrons and muons), assuming total charge neutrality. Thus for a given total baryonic density the proton 
and neutron fractions are determined by the chemical equilibrium condition 

Ain - Aip = 4:l3Esy„i He (10) 

where fip and fj,e are the chemical potential of neutrons, protons and electrons (no muons for simplicity), respec- 
tively. The electron chemical potential is determined by the charge neutrality with protons, assuming the electrons 
to form a free Fermi gas. The crucial property is the density dependence of the symmetry energy which determines 
the imbalance between proton and neutron fractions. In general the nuclear system in such state is strongly isospin 
asymmetric. As a consequence, the calculations of the transport coefficients must be extended to asymmetric nuclear 
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FIG. 3: (color online), (a) Neutron-neutron differential cross sections in pure neutron matter, (b) Neutron-neutron differential 
cross section in pure neutron matter, symmetric and /3-stable nuclear matter, (c) Neutron-proton differential cross section in 
pure neutron matter, symmetric and /3-stable nuclear matter. The free cross section is also plotted for comparison. 



matter. Eq.® for the neutron-proton cross section is to be replaced by the following one 

^ y |G|.,(^)p (11) 

For small asymmetries the isospin effect on the 'reduced' effective mass is of the order, because m* k, itiq + jSra 
and m* w ttiq — /3m, where ttt-q is the effective mass [20l . [2ll | in symmetric nuclear matter. 



IV. NUMERICAL RESULTS AND COMPARISON WITH OTHER CALCULATIONS 

The shear viscosity r] vs. density was calculated for the three different nuclear matter configurations above consid- 
ered, i.e. neutron matter, symmetric nuclear matter, and /3-stable asymmetric nuclear matter. In Fig|4^ the neutron 
and proton viscosities from free and in-medium cross sections are plotted vs. density in the case of /3-stable nuclear 
matter. At low total density the difference between neutrons and protons is about three orders of magnitude because 
the proton fraction is quite small, but it reduces to only one order of magnitude at high density where the proton 
fraction becomes 30% of the total density. As expected, both viscosities, corresponding to free and in-medium cross 
section, are increasing with density, and the medium enhancement turns out to be quite large indeed. In Fig|4jD the 
neutron viscosity in pure neutron, symmetric nuclear and /3-stable matter, is plotted and the medium effects are also 
emphasized by the comparison with the free case. The medium effect is a strong enhancement of 77 , mainly as a 
consequence of the reduction of the level density, more pronounced for /3-stable matter, less for symmetric nuclear 
matter. Comparing neutron matter and /3-stable matter, the values are very close to each other in the low density 
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range, where /3-stable matter is essentially made of neutrons, then they deviate for the increasing weight of the proton 
fraction, according to the /3-stability condition (see Figl3]and Eq.(lO)). 

FigHh shows the comparison among various microscopic calculations. The data from Benhar et al.Q, available up 
to a maximum range 0.3/to~'^, differ from the present calculation for the different neutron-to-proton composition of 
/3-stable nuclear matter. In that case the proton fraction, at the same total density, is larger, as a consequence of a 
larger symmetry energy (curve CBF in Fig[T]), and thus the influence on the neutron viscosity by the neutron-proton 
cross section turns out to be bigger. The data from Yakovlev et al.[Tl| were obtained from the constant effective 
mass approximation (m* jm = 0.8 in the plot), which is not viable for the effective mass gradual quenching at higher 
density. The result is that the neutron viscosity is underestimated and it turns out to be smaller than the electron 
viscosity. The opposite happens in the present as well as the Benhar et al. calculation, as shown in FigHJ:. Therefore 
the electron viscosity would definitely be immaterial in the study of the energy dissipation. 




p [ fm"' ] p [ fm"' ] P [ fm" ] 



FIG. 4: (color online). Shear viscosity from the Brueckner theory, (a) Neutron and proton shear viscosity in /3-stable nuclear 
matter, (b) Neutron viscosity in all considered states, (c) Comparison of the Brueckner neutron viscosity with other calculations 
(see text) and with the electron viscosity, in /3-stable nuclear matter. 



The thermal conductivity k was also calculated according to Eq.(3). The results in the various nuclear matter 
configurations are reported in Fig|SJ Again the medium corrections are emphasized by plotting the k values obtained 
with in-vacuum and in-medium cross sections. In the right panel the present results are compared with the two above 
mentioned models |3, [ll] • To assess the extent of variations it must be taken into account that in this case the scale is 
quite different from that of Fig|4]:. In any case the good agreement among the different models for the pure neutron 
conductivity confirms that the deviations can only be traced to the different isospin dependence. 

The transport parameters, viscosity and thermal conductivity, calculated for /3-stable nuclear matter can be used 
to determine the respective time scales of energy dissipation Ty and tt ■ This requires to integrate ?7(p) and /v(p) 
weighted with the density profile pir) of a given neutron star configuration. The latter can be obtained solving the 
Tolman-Oppenheimer-Volkov (TOV) equation with a given equation of state of nuclear matter. Simplified expressions 
for the time scales governing the energy dissipation from non radial oscillations can be derived for a quasi-uniform 
density model (24l |: 

where R is the radius of the star, p is the density calculated as the ratio between the mass M and the volume AmB? j'i. 
The parameter / is the angular momentum of the non radial oscillation Yim (below we consider l=m=2). As discussed 
in Ref . (24| . Eqs.(12)-(13) underestimate the values for the transport parameters [111 within a factor about 10, because 
the constant mass approximation smoothes out their increase at higher densities. In Table I the damping time scales 
of Y22 non radial modes are reported for a number of neutron star mass-radius configurations. The neutron star 
configurations are taken from Ref. [28] , where the TOV equation was solved using the equation of state derived from 
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FIG. 5: (color online), (a) Thermal conductivity from the Brueckner theory in all considered states, (b) Comparison of the 
thermal conductivity in pure neutron matter between the BHF approximation and two other recent calculations(see text). 



the same BHF approximation as for the in-medium cross sections. The low temperature results show that, within 
the above mentioned uncertainty of constant mass approximation, the time scale for energy dissipation due to the 
neutron shear viscosity could be of the same order of magnitude as the time scale associated to gravitational waves 
(10 — lOOsec), at least at high density. The proton contribution is much less important, being the proton fraction 
much less than the total density. In other words the instability driven by gravitation radiation could be prevented by 
the shear viscosity dissipation in cold neutron stars. To confirm such a statement, a more accurate calculation is to 
be performed by means of the viscosity integration weighted with the neutron star density profiles p(r),0 < r < R. 
On the contrary, the high temperature time scales are too large by several orders of magnitude (ry w lO^'^ 
at T = 3 • lO^^K) and no help can be expected by superfluidity since the critical temperature is much less that the 
newborn star temperature. 

The energy dissipation due to thermal conductivity tt can be calculated from Eq.(12) for a non radial mode I = 2. 
At temperature T = lO^K it is about lO^ry and it gets increased at higher temperature. Therefore its effect on the 
damping can be neglected in agreement with other calculations [24]. 



V. CONCLUSIONS 



The transport parameters, shear viscosity and thermal conductivity, of neutron stars have been calculated in the 
framework of the BHF approximation with two and three body forces. Both forces are described, in a unified treatment, 
by the one-boson exchange model with the meson parameters of the realistic Bonn B potential. In the Brueckner theory 
the in-medium NN cross sections are calculated replacing the in vacuum scattering amplitude with the G-matrix and 
the nucleon bare mass with the effective mass. The in-medium strong effective mass renormalization, which affects 
the level density in the entrance and exit channels, is mainly responsible for the strong deviation of the in-medium 
NN scattering cross sections from the scattering in free space. The main result is a remarkable enhancement of the 
transport coefficients. This effect is well known since some decades, but only recent ab initio calculations provide 
reliable quantitative estimates in domains of nuclear matter, which do not benefit from direct empirical constraints. 
The underlying many body approaches, in fact, are based on realistic NN interactions without free parameters, and 
are able to calculate on the same footing both the equation of state of nuclear matter determining the composition of 
neutron stars as well as the transport properties parameters determining the neutron star cooling and the damping 
of collective motions. 

The calculation of the transport parameters was first performed for pure neutron matter and symmetric nuclear 
matter, then it was extended to /3-stable nuclear matter for the sake of application to neutron stars. The calculation 
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MM@ R(Km) p(grW)- 10" tv(sec) 



0.8 12.8 

1.4 12.5 

1.8 12.2 

2.3 11.0 



1.8 1865 

3.4 1680 

4.7 883 

8.25 317 



TABLE I: Non-linear mode (l=m=2) damping time scale for energy dissipation due to neutron viscosity in neutron stars. 



covers a wide density range, as requested by the study of the neutron-star core. The numerical resuhs are compared 
with other recent ab initio calculations. Concerning the shear viscosity, the BHF prediction is such that the neutrons 
give a contribution larger than the electrons, in contrast to the Yakovlev et al. calculation where the constant 
effective mass approximation is adopted. On the other hand, the different density dependence of the symmetry energy 
is responsible for the different shear viscosity obtained by Benhar et al.@ in beta-stable configurations. 

A preliminary estimate of time scales of non radial mode damping shows that only in cold stars the dissipation 
times due to viscosity are comparable with the gravitational radiation time, whereas the effect of thermal conduction 
is negligible. But these results should be confirmed by more accurate calculations beyond the constant density 
approximation. Such calculations are on the way. 

In order to interpret the spin and thermal evolution, the neutron stars have been assumed to be in a superfluid state, 
that could deeply influence the their transport properties as well. But the recent calculations of the neutron-neutron 
and proton-proton gaps indicate that the proton ^Sq pairing is present only in a restricted region of low density. 
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where the proton fraction in /3-equihbrium with neutrons is quite small [26|, and the high density neutron gap from 
the channel '^PF2 is so small that it could be easily suppressed by even weak polarization effects |27|]. Despite this 
theoretical uncertainty, the onset of superfluidity is supported by the phenomenology, and its role in the damping of 
collective modes should be clarified. 

For the application of the transport properties to neutron stars, one should also include strange components in the 
/3-stable nuclear matter. In fact the strange particle fractions are increasing with density in the inner core, where they 
can contribute with the 20 per cent to the total baryonic matter f25l|. In the case of hybrid stars the quark matter 
can compete with baryons; therefore the contribution to transport properties of quarks must be considered with a 
consistent treatment of hadron-to-quark transition [28| . The last two research issues, superfluidity and transport 
properties of non nucleonic components are presently under investigation. 
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